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Abstract. 

The Binder cumulant at the phase transition of Ising models on square lattices 
with various ferromagnetic nearest and next -nearest neighbour couplings is determined 
using mainly Monte Carlo techniques. We discuss the possibility to relate the value of 
the critical cumulant in the isotropic, nearest neighbour and in the anisotropic cases 
to each other by means of a scale transformation in rectangular geometry, to pinpoint 
universal and nonuniversal features. 

PACS numbers: 05.50.+q, 75.40.Cx, 05.10.Ln 
1. Introduction 

In the field of phase transitions and critical phenomena, the fourth order cumulant of the 
order parameter £Q, the Binder cumulant U, plays an important role. Among others, 
the cumulant may be used to compute the critical exponent of the correlation length, 
and thence to identify the universality class of the transition. 

The value of the Binder cumulant at the transition temperature, U(T C ), the critical 
Binder cumulant, has received much attention, too, indicating the universality class 
as well |2j. For instance, in the case of the spin-1/2 Ising model with ferromagnetic 
nearest-neighbour couplings the critical cumulant has been determined very accurately 
in numerical work, applying Monte Carlo techniques [3 and transfer-matrix methods 
augmented by finite-size extrapolations to the thermodynamic limit [4j. The resulting 
value, U(T C ) = 0.61069... jl], has been observed to hold also for various other, Ising-type 
models on a square lattice jSlHj- 

Quite recently, it has been demonstrated by Chen and Dohm |H] that the critical 
Binder cumulant may display non-universal features, when the couplings on the square 
lattice are anisotropic. More specifically, they have studied critical effects of the 
anisotropy matrix of the Landau-Ginzburg-Wilson Hamilitonian, consisting of the 
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Figure 1. Sketch of the interactions in the anisotropic Ising model. 



coefficients in front of the leading, second-order gradient term of the Hamiltonian. A 
nondiagonal matrix is shown to imply that the isotropic case cannot be restored by a 
rescaling of lengths, i. e. in choosing suitable aspect ratios on a rectangular geometry. 
Such a scale transformation is possible in the case of a diagonal anisotropy matrix, for 
instance, for nearest-neighbour Ising models with two different couplings along the axes 
of the square lattice |H El E] • 

In this Letter, we shall present results of rather extensive Monte Carlo simulations 
for a variety of anisotropic Ising models on a square lattice, corresponding to diagonal 
and nondiagonal anisotropy matrices. Our data agree qualitatively with the recent 
results by Chen and Dohm |5 a on nonuniversality in the case of a nondiagonal matrix. 
They confirm quantitatively the scale transformation given by Kamieniarz and Blote [3] 
in the case of a diagonal matrix. 

2. Results 

We consider spin-1/2 Ising models on a square lattice with anisotropic nearest and 
next-nearest neighbour interactions. The Hamiltonian reads 



where S X)V = ±1 is the spin at site (x, y), see Fig. 1. In most cases, we take all 
couplings to be positive, i. e. ferromagnetic. In the Monte Carlo simulations, systems 
with KL spins, subject to full periodic boundary conditions, are analysed, where K(L) 
corresponds to the x(|/)-direction. The aspect ratio, r, is defined by r = K/L. 
The critical Binder cumulant is defined by pQ 





taking the thermodynamic limit, L, K — > oo; < M 2 > and < 
and fourth moments of the magnetization, M = ^2, S x ,y/{KL). 
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Figure 2. Critical Binder cumulant U(T C ) for Ising models of size L 2 with J v /Jh = 
4.0, 3.0, and 2.0 (from bottom to top). Open symbols denote Monte Carlo data, full 
symbols, at 1/L = 0, result from the scale transformation 0] as described in the text. 
Lines are guides to the eye. 



We studied especially two types of anisotropy. In the first case, where both diagonal 
next-nearest neighbour couplings vanish, J dl = J d2 = 0, the anisotropy matrix is 
diagonal 0. 

In the second case, where J v = Jh, Jai = Jd, and J d i = 0, one encounters 
nondiagonality 5\. Indeed, we consider the two-dimensional variant of a three- 
dimensional model which had been analysed by Chen and Dohm |H] and subsequently 
investigated using Monte Carlo techniques [7j. The simulations, however, may have 
been hampered by the fact that the critical phase transition temperature is not known 
exactly, in contrast to the present situation. 

Note that data of high accuracy are needed. We computed systems of sizes, K, L, 
ranging typically from 4 to 64 for square lattices, K = L. Using the standard Metropolis 
algorithm (a cluster flip algorithm becomes significantly more efficient for larger system 
sizes), Monte Carlo runs with up to 5 x 10 8 Monte Carlo steps per site, for the largest 
systems, were performed, averaging then over several, about ten, of these runs to obtain 
final estimates and to determine the statistical error bars, shown in the figures. We 
computed not only the cumulant, but also other quantities like energy and specific heat, 
to check the accuracy of our data. Of course, sufficiently small lattices may be solved 
easily by direct enumeration. 



2.1. Diagonal anisotropy matrix 

We consider the anisotropic Ising model Eq. (1), with vanishing next-nearest neighbour 
interactions and different horizontal, Jh, and vertical J v , nearest-neighbour interactions, 
see Fig. 1. Then the anisotropy matrix is diagonal [5 . The exact transition temperature 
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is given by [HI E] 

smh(2 J h /k B T c ) sinh(2 J v /k B T c ) = 1, (3) 

where k B is the Boltzmann constant. Following Kamieniarz and Blote jlj, the critical 
Binder cumulant for the anisotropic model on a square lattice, K — L, in the 
thermodynamic limit, can be obtained from that of the isotropic case, Jh = J v , on 
a rectangular lattice with the aspect ratio r = K/L using a 'scale transformation': 

smh(2J h / k B T c ) = 1/ sinh(2 J v /k B T c ) = r. (4) 

The critical cumulant U (T c ) of the isotropic case with arbitrary aspect ratio r has 
been calculated before, extrapolating to the thermodynamic limit, with a high degree 
of accuracy resulting in a polynomial representation of U(T C ) in r jlj. Note that this 
critical cumulant decreases monotonically from its maximal value at r = 1 to zero, when 
the aspect ratio is increased to infinity [I] (or lowered to zero, because of the obvious 
symmetry when replacing r by 1/r), i. e. when approaching the one-dimensional limit 
of the model. 

In Fig. 2 we display results of our check on that scale transformation, showing 
our Monte Carlo results of U(T C ) for anisotropic Ising systems, Jh ^ Jv, of small and 
moderate sizes, K = L. We also include the values we obtained by evaluating the 
previously reported polynomial results for the isotropic case with the appropriate, see 
Eq. (4), aspect ratio r |4.. We find a very good, quantitative agreement between our 
Monte Carlo data and the values based on the scale transformation and the polynomial 
representation jlj. 

This observation confirms the universality of the critical Binder cumulant in the 
case of a diagonal anisotropy matrix. 

2.2. Nondiagonal anisotropy matrix 

Now, we consider the case where J v = Jh = J, and Jd\ = Jd, while Jd2 = 0, see 
Fig. 1. The model has been called, albeit being defined on a square lattice, the 
'anisotropic triangular model' 5], because of an obvious isomorphy. Its anisotropy 
matrix is nondiagonal [Sj. Again, the exact transition temperature is known [9} ITO] 

(sinh(2 J/k B T c )f + 2 sinh(2 J/k B T c ) smh(2J d / k B T c ) = 1, (5) 

where J and Jd are supposed to be ferromagnetic in the following, unless stated 
otherwise (a weakly antiferromagnetic Jd does not destroy the ferromagnetic order at 
low temperatures [Oj). Some of our simulational results for the Binder cumulant at T c 
for square lattices, L = K, with ferromagnetic couplings are shown in Fig. 3, together 
with highly accurate estimates based on finite-size extrapolations to the thermodynamic 
limit |3] in the cases Jd = 0, being the isotropic nearest neighbour square Ising model, 
and Jd = J, being the isotropic nearest neighbour triangular Ising model. The Monte 
Carlo data seem to allow a smooth extrapolation to the thermodynamic limit, providing 
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Figure 3. Critical cumulant of the anisotropic triangular Ising model for sizes L 2 , 
with Jd/J= 0.0 (squares), 1.0 (triangles up), 1.5 (triangles left), and 2.5 (triangles 
down). Open symbols denote our Monte Carlo (as well as numerically exact data for 
small lattices), full symbols, at 1/L = 0, denote previous estimates following from 
finite-size extrapolations 0]. 



reliable and accurate estimates (we refrained from a quantitative finite-size-analysis, 
because, in general, the form of the corrections to scaling seems to be unknown). 

Results for extrapolations to the thermodynamic limit are depicted in Fig. 4. 
Most importantly, the critical Binder cumulant U (T c ) is seen to increase from its value 
in the isotropic, Jd = 0, case when introducing the diagonal interaction Jd > 0. 
U(T C ) seems to reach a maximum at about Jd/J= 1.0, decreasing then, crossing the 
isotropic value at about Jd/ J ~ 1.74 (possibly at a/3, as one may speculate), and 
finally approaching the one-dimensional limit, where U(T C ) tends to go to zero, when 
J becomes indefinitely weak compared to Jd- As has been just emphasized, U(T C ) is 
initially, for < J&j J < 1.74, larger than in the isotropic limit. Now, varying the aspect 
ratio and keeping the rectangular geometry in the isotropic Ising model, this can only 
lead to a decrease in the critical Binder cumulant [3]. Therefore, it is not possible to 
obtain the critical Binder cumulant in the anisotropic case, K = L and < J&j J < 1.74, 
by a scale transformation from the isotropic case, in contrast to the situation discussed 
in the previous subsection. Obviously, this finding is in accordance with the analysis by 
Chen and Dohm [Hj. Our findings do not rule out the possibility that special cases of 
the triangular and square models, choosing suitable aspect ratios and anisotropies, may 
be mapped onto each other, having the same critical cumulants j3]. 

Note that the critical cumulant seems to decrease monotonically when taking a 
weak antiferromagnetic coupling J^, as we find in preliminary simulations. 

Perhaps quite interestingly, the critical Binder cumulant displays a non-monotonic 
behaviour when varying the aspect ratio for the anisotropic triangular model, when Jd/ J 
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Figure 4. Critical Binder cumulant, extrapolated to the thermodynamic limit, see 
Fig. 3, for the anisotropic triangular Ising model as function of the anisotropy Ja/ J ■ 
Full symbols denote estimates taken from Rcf. 4 . The solid line corresponds to the 
isotropic case on a square lattice, Jd= 0. 



is sufficiently large. Starting with the square lattice, r=l, and then, say, enlargening 
(or, equivalently due to symmetry, lowering) the aspect ratio, U(T C ) first becomes larger, 
before finally decreasing again as one approaches the one-dimensional geometry. 

In summary, the Monte Carlo simulations show that for a diagonal anisotropy 
matrix, the critical Binder cumulant in the anisotropic and isotropic cases are connected 
by a scale transformation keeping rectangular geometry. In marked contrast, such a scale 
transformation does not exist, in general, in the anisotropic triangular model described 
by a nondiagonal anisotropy matrix, for reasons which have been discussed recently by 
Chen and Dohm, demonstrating non-universal features in the critical Binder cumulant 

Thence, care is needed in identifying universality classes and in estimating phase 
transition temperatures in using the critical Binder cumulant. 
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